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Abstract. Cross- linked actomyosin bundles retract when severed in vivo by laser 
ablation, or when isolated from the cell and micromanipulated in vitro in the presence 
of ATP. We identify the time scale for contraction as a viscoelastic time t, where 
the viscosity is due to (internal) protein friction. We obtain an estimate of the 
order of magnitude of the contraction time r ~ 10 — 100 s, consistent with available 
experimental data for circumferential microfilament bundles and stress fibers. Our 
results are supported by an exactly solvable, hydrodynamic model of a retracting 
bundle as a cylinder of isotropic, active matter, from which the order of magnitude of 
the active stress is estimated. 



1. Introduction 

Cells organise filamentous actin (F-actin) into complex cytoskeletal structures that play 
a major role in determining cell rheology and cell shape, and in generating the forces 
necessary for motility, cell division, or tissue stability [lj. Depending on the nature of 
the actin-binding proteins that nucleate, help polymerise, and cross-link actin filaments, 
the structure of F-actin networks varies widely [2j |3] , from the near isotropic gels found 
in the cell cortex to the linear bundles that form filopodia or stereocilia jU |5]. Passive 
cross-linkers, such as a-actinin or fascin, contribute to the network elasticity, wheras 
bipolar mini-filaments of non-muscle myosin II motors that pull on neighbouring actin 
filaments of opposite polarity are responsible for contractility. 

Cross-linked actomyosin bundles include the circumferential bundles (CBs) found 
near apical cell junctions of epithelia [6], and the stress fibers (SFs) assembled in the 
cytoskeleton of non-muscle animal cells that exert substantial traction forces on their 
environment |7j. They are robust and stable organelles, that can be extracted from 
the cell while maintaining their shape, structure, and contractility, over durations long 
enough to allow mechanical testing P [9j [ID] . Indeed isolated circumferential bundles 
[9] as well as isolated stress fibers [10] contract in vitro in the presence of adenosine 
triphosphate (ATP). 
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Laser photo-ablation stands out as a powerful and versatile tool, widely used to 
reveal the distribution, orientation, and relative intensity of stress in living cells and 
tissues (see [TT1 [T2] for recent reviews). In order to investigate the cytoskeletal basis 
of the mechanical properties of cells and tissues, the retraction of ablated actomyosin 
bundles has been studied quantitatively by several groups [131 EH EH1 EEl EL7] • These 
experiments give unambiguous evidence that the bundles are under tension in vivo, 
since they retract once severed, and that the retraction is an active process, since 
it is blocked by inhibitors of contractility [IHl [161 EB]- The relaxation dynamics is 
generally fitted by a single exponential, that defines a contraction time, whereas the 
initial velocity at the time of ablation gives access to (relative) measurements of the 
preexisting bundle tension. Interestingly, the retraction that follows laser ablation is 
also relevant in (unperturbed) physiological conditions since stress fibers spontaneously 
rupture in vivo [19J. 

We use dimensional analysis to estimate the contraction time of actomyosin bundles, 
and show that the main source of dissipation during contraction is protein friction 
(section [2]). The mechanical properties of retracting stress fibers have so far been 
interpreted with the help of models [TF1 EEl EH EE] that emphasise the underlying 
sarcomeric structure. One exception is a discrete, numerical model based on the 
tensegrity hypothesis that accounts for much of the observed phenomenology [22J. In 
section [31 we introduce and solve a continuous model of a cross-linked actomyosin bundle, 
valid on hydrodynamic length scales large compared to the mesh size of the F-actin 
network, and where contractility is modeled by a constant active stress term [231 El] • 
In section HI comparison with experimental data yields the order of magnitude of the 
active stress of stress fibers and circumferential bundles. 

2. Contraction time 

The relaxation of the length of a stress fiber severed in vivo is well fitted by an 
exponential function of time, with a contraction time r of the order of 1 to 10 s, 
depending on cell type [HI [TF1 EE], and spatial position within the cell [25J. Laser 
ablation of the apical circumferential bundles of epithelial cells yields a somewhat longer 
contraction time of the order of 10 to 100 s [131 ESI [261 1211 • Retracting bundles 
were first depicted as Kelvin- Voigt bodies [16], using a simple phenomenological model 
consisting in a dashpot and a spring in parallel, with a viscoelastic time equal to the 
ratio of the viscosity coefficient to the elastic modulus. However, the origin of the 
viscosity remains controversial. Internal viscosity dominates according to [IT], whereas 
an external drag force is preferred in [20] , implying a surprisingly large cytosolic viscosity 
of the order of 10 Pa s. In [22], the contraction time closely follows the (microscopic) 
viscoelastic time of components of the fiber, a parameter of the model. 

In section 12.11 we estimate the orders of magnitude of the viscosity coefficient due 
to internal protein friction, of the elastic modulus of the bundle, and deduce the order 
of magnitude of the viscoelastic contraction time due to protein friction. In section 
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\'2.'2\ using well established numerical values of the cytosolic viscosity, we further show 
that external viscous drag is negligible, and would lead to a much shorter viscoelastic 
time, at variance with experiment. We next evaluate the poroelastic time that governs 
the permeation of the cytosol through the F-actin network, and discuss its possible 
relevance to the dynamics of contraction. Finally, we calculate the order of magnitude 
of the elastic modulus and viscosity coefficient of stress fibers and circumferential bundles 
(section [273]) . An order of magnitude is hereafter defined by the integer value, obtained 
by truncation, of the decimal logarithm of the quantity of interest. 

2.1. Viscoelastic time 

We propose that the main source of dissipation in cross-linked actomyosin bundles lies 
in the binding-unbinding dynamics of cross-linkers on actin filaments. This dynamics 
leads to an effective (protein) friction, as introduced first in the context of rubber friction 
[28] . and later for actomyosin dynamics in muscle cells [29] . 

In the linear regime, the protein friction force F p exerted on one filament is 
proportional to its relative velocity v (with respect to neighbour filaments), with a 
friction coefficient ( p : F p = —( p v. Since the energy dissipated after unbinding was 
stored as elastic energy in strained cross-linkers, an estimate of the friction coefficient 
is [28]: 

( P = n x k x r x , (1) 

where n x is the average number of attached cross-linkers per filament, k x is the spring 
constant of the cross-linker, and t x is the typical time for bond rupture. We expect 
that several cross-linker proteins, either passive or active, would contribute additively 
to the friction coefficient: ( p = J2x n x k x r x . To our knowledge, the spring constant 
of a-actinin, arguably the most prevalent passive cross-linker in stress fibers, has not 
been measured. For lack of data, only protein friction due to myosin filaments will be 
taken into account. Using the numerical values n x w 10, k x « 0.1 pN nm -1 for myosin 
filaments (30] [31], and equation (CQ), we obtain the order of magnitude of the microscopic 
friction coefficient ( p ~ 10 -3 — 10~ 2 N s m _1 . 

We expect the dissipative stress a p due to protein friction to scale as a p ~ np( p v /A, 
where np is the number of actin filaments in a section of the bundle, A is the area of 
the section, and v is the velocity. In three dimensions, protein friction translates into 
a viscosity coefficient r] p , with o~ p = rj p ^- « Vpj;^ wnere U and L are the contraction 
velocity and the length of the bundle. Assuming for simplicity a constant velocity 
gradient v z — U4-, the local relative velocity on the scale of a filament is v = U l -j^, where 
lp is the typical length of an actin filament (see figured] (B)). Note that although we 
have assumed that filaments are aligned in the z-direction, we may also consider an 
homogeneous, random orientation of filaments. This would modify the prefactors which 
arise from the average orientation of two filaments. Equating the two expressions for 
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a p , we deduce the order of magnitude of the three-dimensional viscosity coefficient: 

V P ~ n F l -j ( p . (2) 

The stiffness of the bundle is dominated by cross-linkers, such as myosin filaments, 
known to be softer than actin filaments: kx ~ 0.1 pN nm _1 [30| |3T] compared to 
kp fa 10 pN nm" 1 [321 S3]- The elasticity of the bundle is characterised by a network 
structure of actin filaments and cross-linkers with long residual time. The network 
has two ends at z = and z = L connected by various pathways in parallel. Each 
pathway has filaments and cross-linkers connected in series, and the bundle elasticity 
is approximated by that of a dominant pathway (figure [1] (C)). When hard and soft 
springs are connected in series, the force acting on the ensemble of two springs for the 
displacement Az is F ~ Az/(kp 1 + k^ 1 ) where kp is a spring constant of a hard actin 
filament and kx is that of a soft cross-linker. When kp 3> kx, the force is dominated by 
kx- The elastic stress a el is of order <r el = E^- fa E^-, where E and AZ respectively 
denote the Young modulus, and the (macroscopic) displacement of the bundle. Denoting 
by lx the typical length between cross-linkers in the network, approximated by the length 
of a myosin filament, we expect the (microscopic) displacement of a cross-linker to scale 
as Az = AZ Since the number of cross-linkers in a section is nxnp, we find for 
the elastic stress a cl ~ nxnp kxAz/A. Since the two expressions for <r el are equal, the 
order of magnitude of the Young modulus reads: 

E fa n x n F ^ k X - (3) 

Using equations (JH13]), we obtain the order of magnitude of the viscoelastic time 
due to protein friction: 

Vp h fA s 

t= eT x Tx ' (4) 

independent of the extension of the bundle. Indeed, experiments suggest that the 
contraction time is independent of the initial width ao [16], and of the initial length 
Iq [S] of the ablated stress fiber, at odds with the prediction r oc /q of the sarcomere- 
based model introduced in [20J. 

The typical lengths of actin and myosin filaments are respectively lp fa 1 fim 
and lx ~ 0.1 /xm [351 EH]. The association/dissociation rates of myosin light chain 
and a-actinin have been measured in stress fibers in vivo in fluorescence recovery after 
photobleaching assays [37J. The turn-over times are of order rx fa 1 — 10 s, consistent 
with the unbinding times measured on single actomyosin bonds under physiological loads 
[381 139] . We obtain for the order of magnitude of the contraction time r the range 

r fa 10 1 - 10 2 s, (5) 

in accord with experimental observations for contraction in vivo after laser ablation, as 
well as in vitro after extraction from the cell [91 ITU] . 
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Figure 1. (Color Online) (A) Schematic image of the contraction of an actomyosin 
bundle. Permeating cytosol is squeezed out of the network during contraction. One 
end of the bundle is fixed at the origin z = 0. The bundle contracts longitudinally 
from the initial length Lq to the final length Lq(\ + e^) (eoo < 0), and radially from 
the initial radius a Q to the final radius ao + (u^ < 0). The velocities of filaments 
whose center of mass positions are located at z\ and zi are denoted as v(z\) — Uz\/ L 
and v{z2) — UZ2/L, respectively. Cross-linkers contributing to the bundle viscosity 
(resp. elasticity) are drawn in red (resp. blue). Importantly, the turn-over of both 
active motor filaments and passive cross-linking proteins gives rise to protein friction. 
The typical distance between the centers of mass of two filaments is Al « lp as 
schematically shown in (B). The abstract network structure is sketched in (C) where 
many pathways connect the two ends (z = and z = L); each pathway contains 
filaments and cross-linkers in series. 



2.2. Poroelastic time 

One obvious alternative source of dissipation is hydrodynamic friction, due to the 
viscosity of the cytosol. At small scale, the longitudinal friction coefficient due to viscous 
drag on an actin filament, modeled as a cylinder of length lp and diameter dp, reads 
(d = 2iir] cytoso \l F / (\n(l F /dp) + 7//), where r] cytoso \ is the viscosity of the cytosol, and 
7// is a dimensionless number of order 10 _1 [40J. Using lp fa 1 /im, dp ~ 10 nm, and 
well-established values for the cytosolic viscosity r/ cytoso i fa 1CT 3 — 10 _1 Pa s [41], we find 
Cd ~ 10 -9 — 10~ 7 N s m _1 , negligible when compared to ( p . At the scale of the bundle, 
with l B ?a 10 /im, d B fa 0.1 /im, we obtain ( D fa 1CT 8 — 10~ 6 N s m -1 ( p . Viscous 
drag against the cytosol may thus be neglected. 

Further, cytosol is squeezed out of the biopolymer network during the contraction: 
permeation also contributes to the dissipation. The permeation of a solvent through an 
elastic network has been studied in the context of chemical gels (see jl2] for a recent 
review of stress-diffusion coupling), and has led to interesting insights into, e.g., the 
mechanics of the poroelastic cytoskeleton [41] . In a two-component system where fluid 
cytosol (the solvent) permeates an elastic polymer matrix, the total stress is expressed 
as the sum of the elastic stress of the polymer network, and of a pressure term for the 
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liquid solvent: aff = erf- — p 5ij, where 6^ is the identity tensor (see section [3] for a 
detailed calculation). 

Supplementing this constitutive equation with Darcy's law for the permeation of 
the solvent, one obtains a diffusion equation for the displacement field of the gel, with a 
diffusion constant D ~ kE, where k is the permeability coefficient. Since k oc £ 2 /?7cytosoi, 
where £ is the mesh size of the network and ?7 cy tosoi is the cytosolic viscosity, the time 
scale for relaxation over a typical distance a , evaluated as the radius of the bundle, 
reads: 



From equation (J3J) , and A ps a 2 ,, we expect the poroelastic time to be sensitive to 
the width of the bundle: The typical radius of a stress fiber is Oq F 100 

nm |43j, whereas the cross-sectional area of a circumferential bundle is of the order of 
A CB n 1 /im 2 El Q31 E], or a£ B « 1 fim. Using £ « 10 nm [TO] and n F « 10 [TO], 



Permeation of the cytosol through the polymer network occurs on diffusive times 
that are shorter than the observed contraction time, and may be safely ignored. It 
may however become relevant for very thick bundles: note that the upper bound of the 
estimate of the poroelastic time is only one order of magnitude shorter than observed 
contraction times in the case of circumferential bundles. In section I3T21 we introduce and 
solve a theoretical model of an actomyosin bundle as an active, viscoporoelastic material. 
We show in section 13.41 that the viscoelastic behavior, as discussed in section I2.1[ is 
recovered in the limit where the viscoelastic time is much larger than poroelastic diffusion 
times. 

2.3. Material properties 

The Young modulus may be directly estimated using equation (jSJ). For stress fibers, 
A SF « 10~ 2 /im 2 , whence E SF ps 10 5 Pa. In (131 H4l . the force-extension curve of an 
isolated stress fiber has been measured in the absence of ATP, up to large deformations 
of the order of 200% where the fiber breaks. In the linear regime of small strains, the 
fiber's Young modulus was estimated to be of the order of 10 5 to 10 6 Pa. Circumferential 
bundles are thicker A CB ps 1 /im 2 , and correspondingly softer: we obtain from equation 
(j3]) E CB p^ 10 3 Pa, somewhat higher than the value E CB p^ 10 2 Pa given in JJ5], where 
the elastic modulus of the actin cortex was measured in vivo in confluent endothelial 
cells. 

Estimates of the bundle viscosity are obtained from equations ([lj and (T5]). Since: 




(6) 



we obtain: 



r p SF oc 10~ 6 - 10~ 4 s, 
r p CB oc 10~ 2 - 10° s. 



(7) 
(8) 



V P ~ n x n F -r k x r x , 



(9) 
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we find: 

r]f ps 10 6 - 10 7 Pa s, (10) 
r£ B ^ 10 4 - 10 5 Pa s. (11) 

Although the effect of internal friction has been considered in models of stress fiber 
retraction [T71 [20], this constitutes, to our knowledge, the first calculation of the 
numerical value of the viscosity coefficient. We emphasise that this value is an order of 
magnitude, derived from dimensional analysis, and defer a more rigorous treatment to 
future work. 

3. Models of active contraction 

We propose a three-dimensional, continuous, exactly solvable model of a cross-linked 
actomyosin bundle as an active, cylindrical body, inspired by active gel descriptions 
of the cytoskeleton [231 EI], where chemo-mechanical coupling is taken into account 
at linear order by including an active stress term in the constitutive equations. The 
resulting dynamical equations can be solved analytically, and yield results that agree well 
with experimental data. Related approaches successfully accounted for the possibility 
of oscillations in muscle sarcomeres, described as active elastic bodies [HH H7] and for 
the existence of distinct types of polarity patterns in contractile actomyosin bundles 
described as active, polar elastomers [38]. The assumptions made are discussed in detail 
in section 13.11 We introduce in section 13.21 a theoretical description of an actomyosin 
bundle as an active, viscoporoelastic material. We discuss in section 13731 the poroelastic 
limit. In section [3741 we show rigorously that the bundle behaves as an active viscoelastic 
body when the viscoelastic time is much larger than the poroelastic diffusion times, and 
solve the model in this limit. 

3.1. General formulation of the model 

We formulate a hydrodynamic description valid on length scales large compared to 
the typical mesh size of the underlying polymer network d 3> £ ~ 10 nm [HI [10]. In 
experiments, bundles are ablated far from their end points. We consider the retraction 
of one half of the severed bundle: the other tip is fixed and sets the origin of the z axis. 
We consider the simplified geometry where (half-) a bundle initially adopts the shape 
of a cylinder of axis Oz (see figure [H (A)). We naturally use the cylindrical coordinates 
(r,9,z). The bundle initially occupies the volume defined by < r < a , < z < Iq, 
where a and l respectively denote the initial radius and length. 

While this simplified geometry may seem a natural approximation for stress fibers 
spanning the ventral side of adherent cells and bound to the substrate at (fixed) focal 
adhesions, we believe that the following observations also justify its relevance for apical 
circumferential bundles. First, the ablation is generally performed along nearly straight 
cell junctions: the local curvature may be neglected. Second, the shape of the cell is most 
often preserved despite the ablation of one cell junction, while other apical bundles in the 



Contraction of cross-linked actomyosin bundles 



8 



same cell remain mostly unperturbed [131 EE] , perhaps due to the stabilising influence 
of the surrounding epithelium. The retracting cell junction may therefore be treated as 
attached to an immobile cell-cell vertex and nearly independent of the remainder of the 
apical actin cortex. 

The severed stress fibers do not retract when the ablation is performed in the 
presence of inhibitors of contractility [101 [16] . For this reason, the displacement field u 
is defined with respect to the initial state (immediately after ablation): u(t = 0) = 0. 
Following [19] , we postulate the following functional form for components of u: 



This Ansatz is central to our calculation. The displacement field is independent of 9 
due to rotational invariance about the axis Oz. The origin is fixed, u z (z = 0) = at 
all times, and the longitudinal strain e(t) is uniform. The initial cylindrical shape is 
preserved by our Ansatz for the displacement field, and the radius a(t) and length L(t) 
of the cylinder at a later time t > read: 



where e(t) appears as the relative variation of length of the bundle, negative for a 
contraction (see figure [1] (A)). 

For simplicity, we assume that the bundle material properties are isotropic. In 
section [21 our estimates rely on the physical properties of cross-linkers, which remain 
approximately isotropic. In addition, the orientation of actin filaments inside actomyosin 
bundles shows a significant amount of disorder. The binding/unbinding kinetics of cross- 
linkers is taken into account by an effective, bulk viscous term with coefficient r) p , due 
to protein friction. We include in the constitutive equations a constant, isotropic, active 
stress <r* ctlve = cr^Sij, and a a is positive since the bundle is under tension. The active 
stress o~a is assumed to be uniform, and independent of time: we neglect the banded 
myosin patterns that reflect the (disordered) sarcomeric organisation typical of stress 
fibers [50] • Note however that not all contractile actomyosin bundles exhibit a sarcomeric 
organisation [5TI 152] . 

On the time scale involved, the final rest state is well defined. For stress fibers, 
the relative total retraction length |eoo| is of the order of 10% [El [IT]: deformations are 
small enough for linear elasticity to hold. We denote by E and v the Young modulus 
and the Poisson ratio of the bundle, and by K = „ / , £ n ; and G = n/ , E , ? the bulk and 

' J 3(1— 2v) 2(l+i/) 

shear elastic moduli. Shear strains and shear elastic stresses cancel by construction. 
The normal elastic stresses read: 



u r (r,9,z,t) 
u e (r,6,z,t) 
u z (r,9,z,t) 



u(r, t) 


e(t) z. 



(12) 
(13) 
(14) 



a(t) = a + u(a , t), 
L(t) = L (l + e(t)), 



(15) 
(16) 




(17) 



(18) 
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Like any organelle in a living cell, actomyosin bundles are highly dynamic 
structures: protein components constantly renew while the global organisation of the 
bundle is preserved. For stress fibers, fluorescence recovery assays performed in vivo 
yield an association/dissociation time of the order of a few minutes for actin j37J 153] - 
In this model, we assume that a steady-state is reached over the duration of ablation 
experiments, and accordingly that the material content of the bundle is constant during 
relaxation. 



3.2. Active viscoporoelastic contraction 



In this section, we study the permeation of a viscous liquid (the cytosol) through an 
active, elastic polymer network with transient cross-links. Active permeating gels have 



recently been studied in detail [54] 155] . Here, we follow a simpler approach, inspired by 
an exact calculation pertaining to the free swelling of a chemical gel whose initial shape 
is a cylinder [19], including two additional ingredients: a constant active stress, and a 
bulk viscous stress due to protein friction. 

Since the strain field obeys the Ansatz (I12til4p . shear strains and shear elastic 
stresses cancel. The normal stresses read: 



0~ r 



^ P ~3r A ^r, 

V r . el 

r 

9v z , el 

Vp ~dz A ° zz ' 



(20) 
(21) 

(22) 



where v = ^ denotes the network velocity field. For a constant active stress a a, force 
balance reads: 



d r p 



d 6 p 
d z p 



d r af r + 



cl 



a rr a 6 1 



3/l_9 
dr \r dr 



(23) 



(24) 
(25) 

where p is the pressure field of the cytosol. Since dgaf e = d z af z = 0, we have 
dgp = d z p = 0: the pressure field is a function of radius and time only, p = p(r, t). The 
radial dependence of the pressure field is obtained from equation (123]) after integration 
with respect to r, using the boundary condition p(a,t) = p ext , where p ext denotes the 
hydrodynamic pressure in the surrounding cytosol. We find: 



p{r, t) = p cxt + r] p 



1 d_ 

r dr 



rv r 



'1 - v)E 



;i + 2i/) 



1 ■ 

- — [ru) 
r or 



(26) 



Using the previous equation, we eliminate the pressure field from Darcy's law V • 
v = k V 2 p, where k is the permeability, and obtain the following equation for the 
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d it + \% r = D ^r^r M ) + S £ {iWr (nv) ) ' (27) 
where Z| = kt^, D = n^rjz^) i s a diffusion coefficient, and their ratio 

Jl = (1 + „)(1 -&,) ,-, 

D l-i/ £ v ; 

is a viscoelastic time. 

The boundary condition at r = a, <7 rr (r = a, t) = — p ext , yields: 

/. ^du ( u(a) \\ dv r , , 

+ H — + e +*ft.-sf. =-^. (29) 



(l + i/)(l-2i/) V dr 

la \ (3 / / dr I a 

In the case of free contraction, the external force applied on the tip of the cylinder 
at z = L is equal to Jq dr l-nr o~ zz (r,t) = —na 2 p ext . Combining the two boundary 
conditions, we obtain: 

d / u(a) \ / u(a) \ , . 

Tg dt V <i ~ J ~ ~ \<* ~ } ' <30) 
where tq is a viscoelastic time based on the shear elastic modulus: 

ra = (31) 

which contributes to the relaxation provided that "( a °'*~°) ^ e (t = 0), see also equations 
(1A.12HA.3~3l) . Here, the initial strain is equal to zero u(a ,t = 0) = e(t = 0) = 0, and 
the solution of equation ( 1301) is for all time t: 

e(t) = ^ UM . (32) 
a ao 

The length and the radius of the bundle both relax in time according to the same 
functional form since e(t) = L<yt \~ L ° = a W~ a ° . 

Equations (I27ti291l32|) are solved by the expansion: 



e(t) -eoc+Eeue"^, (33) 



71 



u(r,t) = Uoo(r) + ~Tu n (r) e" i/T ", (34) 



where the r n are relaxation times. Equation (jzTj) yields a differential equation obeyed 
by the amplitudes w n and e n : 

u n + \e n r = D (^i) ±(l± (ruSj • (35) 

In the following, we define the lengths l n by l\ = D \r E — r n \. 
In the long time limit, the stationary solution reads: 

eoo = - (36) 
«oo(r) = - ^ r. (37) 
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Remarkably, we predict that actomyosin bundles also contract along the radial direction 
(doo < do), independently of the value of the Poisson ratio u, since the final radius reads: 

= a (1 + e^) = a - . (38) 

Indeed, electron microscopy imaging indicates that isolated stress fibers contract radially 
[TO] . To our knowledge, this observation has not been confirmed by optical imaging in 
a live cell, due to insufficient spatial resolution. 

The relaxational dynamics are solved as follows. When te > r n , the functions u n (r) 
are the solutions of an inhomogeneous, modified Bessel equation: 

u n (r) = b n I x (J-^ - ^e n r, (39) 

where the coefficient b n is an amplitude, and h(x) is a modified Bessel function of the 
first kind. Since u n (ao) = a§e n (see equation (1321). we obtain e n = |^ I\ {^) from 
equation fl39|) . We deduce from equation (1291) that the relaxation times r n are solutions 
of the equation: 



^-1 



^ Pl (p)-l h (p. 



J 1 V / / S \ / 



4 " h(?). («) 



3 1 - v \l 



When r n > r E , equation ( )35|) is an inhomogeneous Bessel equation. Modified Bessel 
functions are replaced by Bessel functions: 

u n (r) = b n Ji (j^J - ^e n r, (41) 



where J\{x) is a Bessel function of the first kind, and e n = |^ J\ yf 1 )- The relaxation 
times T n are solutions of the equation: 



^-1 



4 " (42) 



3 1 - i/ V/ 



/ 1 \ 1 3 W 

which may be obtained from equation (T4Tjl) upon replacing Ji by J\. In figure EJ we plot as 
a function of the parameter — the dimensionless relaxation times — , at a given value 
of the Poisson ratio. Equation ( 14T)1) admits only one solution, whereas equation ( 1421 ) 
admits an infinity of solutions, which contribute to the expansions ( 1331I34"]) . 

3.3. Active poroelastic limit 



In the limit where r n 3> re, equation (1421) becomes: 

fl ° 4 f-1 + Ji (-) = 0, (43) 



r n \r n J 3(l-i/) \r n 
with r ra = y/Dr n , in agreement with jl9]. The largest solution of equation ( f43|) is 

a 2 

r n K, t p = this limit is realised when r p 3> te, as may be possible, e.g., for very thick 
bundles assembled in vitro, since we expect the poroelastic time to scale as the fourth 
power of bundle radius, see equation (JBJ) and below. 




Figure 2. (Color Online) Dimensionless relaxation times — as a function of the 
parameter — , obtained numerically for v — 0.4. The dashed blue line (— = — ) 

Tp Tp Tp 

separates the two sets of solutions. The black line gives the unique solution of 
equation (|40l) . which converges to ^ as — > oo (active viscoelastic limit, dashed 
green line). The red lines are the solutions of equation (j42|) . and converge to constant 
values as — — > (active poroelastic limit, solutions of equation (I43[) V 



3.4- Active viscoelastic limit 

We formally define the viscoelastic limit as the regime where the permeability coefficient 
k diverges. Darcy's law then implies that the cytosolic pressure field is homogeneous 
(Vp = 0). Since ao/l n oc we ma y use xl[(x) ~ I\{x) and xJ[(x) ~ Ji(x) in 

equations ( HOI) and ( 142]) since < i < 1. For both equations, we obtain a unique 
solution r n — )■ tk- 

r K = A, (44, 



a viscoelastic time based on the bulk modulus K. The expansions (I33H34I) reduce to one 
exponential term. Given the initial conditions u(r, t = 0) = e(t = 0) = 0, we find the 
expression of the displacement fields: 



«rM) = -||r (l-e-*^), (45) 

u-CM) = (l-e-*/^). (46) 

The length of the ablated bundle decreases as an exponential function of time, as 
observed in experiments. Equation ( H6l) also yields the initial retraction velocity Vq, 
proportional to the active stress: 

v = v z {z = L , t = 0) = -L — . (47) 

?7p 



We may rewrite the argument ao/l n as a function of time scales: p = Jyf^z^r-y 
We find that — = jz^Tk'- clq «C Z n corresponds to r n — > r A - » r s — r A - ^> r p . The 
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contraction of an active viscoporoelastic cylinder (with protein friction) does indeed 
reduce to an active viscoelastic behavior, provided that the viscoelastic time is large 
compared to the poroelastic time. The cross-over between the two regimes occurs when 
T P ~ t k, or <2q ps kt] p : it may be controlled by the bundle radius, with a cross-over 
radius a c . Using k ps £ 2 /?7 cy tosoi and r\ v ps n F ^ ( p (equation ([2])), we find: 

ac ^( n F i F c P e \ ll \ (48) 

\ ?7cytosol / 

Using the same numbers as in section [21 we obtain the value a c ~ 10 /mi, large compared 
to the radius of actomyosin bundles found in cells. The separation of time scale r K 3> r p 
no longer exists when a 3> a c , and is absent without protein friction, i.e. when cross- 
linkers are fixed and do not unbind over the time scale considered. 



In |Appendix A[ we solve a model of an actomyosin bundle as an active, viscoelastic 
material, and check that the results are identical to the limit presented in this section. In 
addition, we consider the case where the initial state differs from the reference state, and 
find that the strain, initially nonzero, relaxes as a linear combination of two exponentials 
with characteristic times t k and r G , equations (1A.12HA.13j) . 



4. Active stress 

In experiments, the stationary longitudinal strain is of order je^l ps 10 _1 . Assuming a 
value of the Poisson ratio of order v = 0.4 [56j [57], the estimates of section [2] remain 
relevant since 3K ps E. Using equation ( |36l) . we deduce the order of magnitude of the 
active stress a a = 3K\eoo\ ps -E|eoo|: 

of ps 10 4 Pa, (49) 
ps 10 2 Pa. (50) 

Since the stall force of a myosin motor is of order F$ ps 1 pN [58] , the above values agree 
with a naive estimate of the active stress as 

o A ps n x n F — , (51) 

where n\ is the number of active myosin molecules (per filament) that generate stress, 
with the same order of magnitude as the total number of attached cross-linkers per 
filament ~ fix. Since L — = |eoo|L , combining equations (JH]) and (|5Tjl yields: 

F s 

Lq — ~ ~, — r~ Lo- (52) 
k x lx 

Experiments suggest that the total retracted length is independent of the bundle radius 
a US]. 

The value of the active stress obtained for circumferential bundles is close to the 
order of magnitude obtained for cortical actin o^ A ps 10 3 Pa [59]. Estimate ( 14"9]) suggests 
that cells whose stress fibers pull on an area of the order of 1 fim 2 exert locally a force of 
the order Fa ~ 10 nN. This is consistent with measurements of traction forces exerted 
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by cells on micropatterned pillars [60J . Adherent cells assemble thicker and more robust 
stress fibers when the rigidity of the substrate is large [61J . In this case traction forces 
saturate to a value independent of the external stiffness. In |62j, a simple model of the 
cytoskeleton as an active, elastic material allowed to identify active contractility as the 
physical origin of the saturation traction force, yielding 10 4 Pa as an order of magnitude 
of the active stress, in agreement with ( I49p . 

The values of measured in ablation assays are similar to the values of the 
"preexisting strain" measured in isolated stress fibers [43 J, or of the "preextension" upon 
unloading the stress fibers of cells adhering on a stretchable substrate [63]. Indeed, 
as we have seen in fl36H37p . actomyosin contraction may be seen as equivalent to a 
reduction in the stress-free reference length of an elastic material [64J. Experiments 
show that the amplitude of the preextension correlates positively with contractility 
[63] , in qualitative agreement with equations ( l36ll3~Tl) . Further, knock-down of a-actinin 
increases preextension [63J. Given that stress fibers depleted in a-actinin are expected 
to be softer [65], this observation also fits with equations ( 13611371) . since a smaller value 
of the elastic modulus at fixed active stress leads to a larger strain. 

5. Conclusion 

Local ablation is a widely used tool to estimate the value of the local stress in living 
cells, up to a (generally unknown) viscosity coefficient (see equations ( ]4"T1) and ( 1A.14I) ). 
Ascribing the physical origin of the fiber's viscosity to protein friction, itself due to the 
association/dissociation dynamics of cross-linkers, we obtain the order of magnitude 
of the viscosity coefficient r) p . Using an independent estimate of the bundle's elastic 
modulus, we deduce the order of magnitude of the contraction time r « 10 — 100 s, 
interpreted as a viscoelastic time. We find that the order of magnitude of the poroelastic 
time scales is much smaller, and study rigorously the limit where poroelasticity may be 
neglected thanks to an exactly solvable model of the contracting bundle. We identify a 
well-defined threshold, controlled for instance by the bundle width: viscoelasticity (resp. 
poroelasticity) dominates when the viscoelastic time is much longer (resp. much shorter) 
than the poroelastic time. For an isotropic material, the model predicts that contraction 
occurs both in the longitudinal and radial directions, irrespective of the value of the 
Poisson ratio. From the observed longitudinal strain, we deduce the order of magnitude 
of the active stress: <r^ B ~ 10 2 Pa in circumferential bundles and <r^ F ~ 10 4 Pa in 
stress fibers. Following the literature, we distinguish stress fibers from circumferential 
bundles. However their protein constituents and macroscopic properties are identical: 
the differences we emphasise may in fact be superficial, and turn out to conceal a 
continuum of parameter values. 

We expect our results to be relevant to the mechanics of other cytoskeletal 
structures, such as the actin purse string at the circumference of healing wounds [66] , or 
the supra-cellular actin cable that surrounds the amnioserosa during dorsal closure of 
fly embryos [67]. Circular, contractile actomyosin bundles also form in adherent (single) 
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fibroblasts during the course of spreading (SB]- Modeling a retracting bundle as a straight 
cylinder is relevant locally along a ring as long as the retracted length is small compared 
to the radius. It may be interesting to generalise the calculation presented here from 
a cylindrical to a toroidal geometry. Strikingly, ablation experiments performed at the 
tissue level in Drosophila pupal epithelia, on the scale of one hundred cells, reveal a 
relaxation time of the same order as the contraction time of individual circumferential 
bundles [69], suggesting that the same physical mechanism may be at play. 

We strongly believe that estimates of orders of magnitude, based on dimensional 
analysis, are a useful contribution to quantitative biology [70]. More quantitative data 
will be needed in order to reduce the spread of numerical values. One promising 
direction in which our ideas may be tested is that of bundles reconstituted in vitro 
[7T| l72| l73l [74] . where microscopic and macroscopic parameters may be varied in a 
controlled fashion. For instance, we expect that the contraction time increases linearly 
with the average length of actin filaments lp (equation (jl])), a prediction that may be 
tested experimentally. On the theoretical side, one would like to be able to calculate 
in a rigorous manner the material parameters of cross-linked actomyosin bundles from 
the microscopic parameters pertaining to individual constituents (filaments and cross- 
linkers). This important question is left for future study. 
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Appendix A. Active viscoelastic contraction 

Using the same Ansatz for the displacement field, equations ( fT2lfl4l) . the constitutive 
equations for the active, viscoelastic contraction of a cylinder of isotropic material read: 

Ov E ( Ou Vj \ 

arr = Vp _ r + aA+ {l + v){l _ 2v) (a--)^+ K" + e)j, (Al) 

v r E ( . .u , du . \ , . „. 

r (1 + v)(l — Iv) \ r or J 

de E ( , . ,u du\ , . „. 

g " = "-di + ^ + (l + ,)(l-2„) + (A.3) 

Again, since u and e are independent of 9 and z, the azimuthal and longitudinal 
components of the force balance equation are immediately verified. In the radial 
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V P d r (rv r )) + {l ^~ V } 2v) dr (h (ru)) = 0. (A.4) 
A radial displacement linear in r: 

u(r, t) = -A(t)r (A.5) 

solves equation flA.4[) . provided that the dimensionless function A(t) depends on time 
only. The stress field is then uniform (with a ee = a rr ): 

1 dA E /.„ .A ,A \ 

*rr = ~ 2 % ^ + O-A + {l + u){l _ 2v) " ") 2 + "( 2 + C) J ' (A ' 6) 

a - = ^l + ^ + (i + ,Ki- 2 ,) « 1 ~ 1/ ) e+ ( A - ? ) 

and determined by the boundary conditions. 

Neglecting the external, cytosolic pressure, normal stresses at the free surface are 
a nn \ r=a = a nn \ z= l = 0. We obtain a system of two coupled differential equations for 
A(t) and e(t): 

(1 + - 2i/)| ^ = - A - 2ve - 2(1 + v){l - 2v)^, (A.8) 

(1 + u)(l - 21/) | f f = - vA + {y - l)e + (1 + - 2z/)^, (A.9) 
diagonalised as: 

r^(A + e) = -(A + e)-££, (A.10) 

r G A(2e-A)= -(2e-A), (All) 

where tq and t# are defined in equations fl3T|) and (|44p respectively. For zero initial 
displacements, we find A(t) = 2e(t), and e(t) = — |^ (l — e~*/ TK ), in agreement with 
equations fl45M6p . 

If we now take into account a possible difference between the initial state and 
the stress-free reference state of the passive bundle, the initial conditions become 
A(t = 0) = Aq 7^ 0, e(t = 0) = eo 7^ 0. The strain relaxes as a linear combination 
of two exponentials with characteristic times tq and tk'- 

A(t) = -^(l-e-^«) + 2 -(A + e )e-^ 

+ ±(A -2e )e- t / Ta , (A. 12) 



e[t) = ~^V~ e ) + 3 ^ + e °) e 

-i (A -2e )e- t ^. (A.13) 
A nonzero initial elastic stress also contributes to the recoil velocity: 
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with o~f z (t = 0) = ( 1+ „)^ 1 _ 2l/ ) ((1 ~~ z/ ) e o + ^o) • Since ^ = the active contraction 
of bundles initially stretched or compressed may yield a measurement of the Poisson 
ratio. In [H], the relaxation of ablated circumferential bundles has been fitted by a sum 
of two exponential functions with relaxation times T\ and r 2 , with ^ « 10. Assuming 
that our model is valid, this observation suggests that, in this experiment, ^ ~ 10, or 
equivalently ^ — v « 10 _1 . 

For completeness, we also solve the case of an incompressible material. The 
dynamical equations for Ait) and e(t) read: 

+ (A.15) 

We find that the bundle contracts both radially and longitudinally with a characteristic 
time To-, the viscoelastic time based on the shear elastic modulus: 

Mt)=-%+C% + Ao)e-'"°, (A. 17) 
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